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ABSTRACT 

It has been suggested that the observed rotation periods of radio pulsars might be induced by a 
non-axisymmetric spiral-mode instability in the turbulent region behind the stalled supernova bounce 
shock, even if the progenitor core was not initially rotating. In this paper, using the three-dimensional 
AMR code CASTRO with a realistic progenitor and equation of state and a simple neutrino heating 
and cooling scheme, we present a numerical study of the evolution in 3D of the rotational profile 
of a supernova core from collapse, through bounce and shock stagnation, to delayed explosion. By 
the end of our simulation (~420 ms after core bounce), we do not witness significant spin up of the 
proto-neutron star core left behind. However, we do see the development before explosion of strong 
differential rotation in the turbulent gain region between the core and stalled shock. Shells in this 
region acquire high spin rates that reach ~150 Hz, but this region contains too little mass and angular 
momentum to translate, even if left behind, into rapid rotation for the full neutron star. We find also 
that much of the induced angular momentum is likely to be ejected in the explosion, and moreover 
that even if the optimal amount of induced angular momentum is retained in the core, the resulting 
spin period is likely to be quite modest. Nevertheless, induced periods of seconds are possible. 
Subject headings: Hydrodynamics — supernovae: general — pulsars: spin 



1. INTRODUCTION 

One of the key questions associated with radio pulsars 
is the origin of their spin periods (Lorimer 2009,2010). 
In two recent papers (Blondin & Mezzacappa 2007 
[hereafter BM07]; Blondin & Shaw 2007), it has been 
suggested that the growth in 3D of an £ = 1 non- 
axisymmetric (m = { — 1,1}) spiral mode (Fernandez 
2010) of the so-called standing accretion shock instabil- 
ity ("SASI") during the post-bounce delay phase seen 
in current supernova theory can spin up matter in the 
turbulent region behind the stalled shock. BM07 and 
Blondin & Shaw (2007) find that the accretion of this 
spinning matter (and its associated angular momentum) 
onto the inner core can spin up the nascent neutron star 
to rotation periods that compare favorably with the ob- 
served/inferred values for pulsars at birth. More pre- 
cisely, they report final periods of ~70 milliseconds (ms) 
for moderately rotating progenitors and ~50 ms for non- 
rotating progenitors. The measured period range of birth 
periods is from ~10 milliseconds to seconds, with an av- 
erage "injection period" near ~0.5 seconds (Chevalier & 
Emmering 1986; Narayan 1987) 3 . If this mechanism were 
to obtain, it would obviate the need to invoke a spinning 
progenitor core. However, rotation is a natural feature 
of stars, and there is no obvious reason the progenitor 
"Chandrasekhar" mass in the center of a massive star 
at death wouldn't be rotating with the spins and an- 
gular momenta needed to leave a compact object with 
the measured spin frequencies (Ott et al. 2006ab). This 
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is a function of the angular momentum distribution in 
the cores of massive stars when they collapse (Heger, 
Woosley, & Spruit 2005; Maeder & Meynet 2000,2004; 
Hirschi, Meynet, & Maeder 2004). Hence, the possibility 
that the 3D dynamics of the core itself could generate 
countervailing flows with opposite signs of angular mo- 
mentum, "one stream" of which could leave the residue 
spinning at interesting rates while conserving overall an- 
gular momentum, is intriguing. If this were true, to ze- 
roth order all cores at collapse could be non-rotating, and 
the spins could be a function of collapse dynamics alone. 
This would eliminate one of the current ambiguities in 
massive star progenitor models, whose core spins at col- 
lapse are not observable. The progenitor cores could all 
be born with "zero" angular momentum, while giving 
birth to proto-neutron stars with respectable spins. 

However, BM07 and Blondin & Shaw (2007) employed 
various simplifications in their simulations. They did not 
follow collapse itself, but started in a steady-state post- 
bounce configuration. They used a simple "gamma-law" 
equation of state (EOS), and not a full nuclear EOS. 
They set the mass accretion rate through the shock to 
a constant value and did not perform their simulations 
in the context of a realistic massive-star progenitor. Im- 
portantly, they excised the inner core whose very spin-up 
was being studied. Rather, they inferred from ostensible 
accretion of matter and angular momentum at their in- 
ner boundary an irreversible accumulation of spin in the 
proto-neutron star. Finally, they did not consider the 
potential role of the supernova explosion itself on the 
magnitude of the residual angular momentum and final 
induced core period. They did, however, investigate the 
potential role of initial rotation in the final outcome. 

In this paper, we investigate the claims of BM07 and 
Blondin & Shaw (2007), improving upon their study in 
several ways. Our 3D simulation is described in Nord- 
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haus et al. (2010a) and was carried out using CASTRO 
(Almgren et al. 2010), an Eulerian Adaptive Mesh Re- 
finement (AMR) hydrodynamics code with hierarchically 
nested rectangular grids that are refined both in space 
and time. We employed a realistic EOS (Shen et al. 
1998ab) and a simple neutrino transfer scheme (see Nord- 
haus et al. 2010a and Murphy & Burrows 2008). Both 
Nordhaus et al. (2010a) and Murphy & Burrows (2008) 
incorporated the Liebendorfer (2005) prescription for 
electron capture on infall and a simple neutrino heating 
and cooling algorithm after bounce. The driving electron 
and anti-electron neutrino luminosities were kept con- 
stant at 1.9 x 10 52 ergs s -1 . This enabled us to heat the 
so-called "gain region" (Bethe & Wilson 1985) and inves- 
tigate the pre-bounce infall, bounce, post-bounce delay, 
and explosion phases. In this simulation (model 3d:L_1.9 
of Nordhaus et al. 2010a), a neutrino-driven explosion 
occurred approximately ~200— 250 ms after bounce and 
we carried the simulation out to ~420 ms after bounce. 
Our initial model was the 15-M solar-met allicity, non- 
rotatin g red-supergiant model of Woosle y" fc Weaver | 
(1995). An important aspect of our simulation is that 
the proto-neutron star itself is at all times present on 
the grid (which is Cartesian) and is followed throughout 
the run in a hydrodynamically consistent fashion. 

The cubic grid in our simulations extended to 10 4 kilo- 
meters (km) on each side, with 304 3 cubic cells covering 
the grid's volume at its coarsest level. Three levels of re- 
finement of a factor of 4 in each dimension were followed. 
As a result, the resolution in the central regions (the in- 
ner 200 km) was ~0.5 km. The monopole approximation 
for self-gravity was adopted. 

By the end of our simulation, we find that there is 
almost no spin up of the proto-neutron star. During 
the post-bounce phase and before explosion, the region 
exterior to the core and interior to the shock (~60 to 
~250 km) does witness the growth of spiral modes, some 
spin up (to ~150 Hz), and the emergence of countervail- 
ing rotational flows. However, the explosion ejects not 
only mass, but its associated angular momentum, and 
depending upon where the mass cut is finally imposed, 
the residual core can end up with some net rotation. In 
our most optimistic scenario the core achieves a period of 
only 1.2 seconds. This fastest possible spin does not alone 
seem adequate to explain the measured and/or inferred 
pulsar birth periods without significant initial rotation in 



the progenitor itself. In { 



we present our quantitative 



results and we follow in [3 with an extended discussion 
and our summary conclusions. 

2. RESULTS 

2.1. Profiles and Temporal Evolution of Angular 
Momentum 

We begin our analysis by calculating the amount of an- 
gular momentum present in our computational domain. 
We divide the innermost 2000 km of the grid into con- 
centric spheres of increasing radius, and we monitor the 
magnitude of each of the three components (x, y, and z) 
of the angular momentum (Li) enclosed within each 
sphere for the duration of the calculation after core 
bounce. The panels in Fig. [I] depict three-dimensional 
surface representations of the temporal evolution of the 
three orthogonal components of the total angular mo- 



mentum enclosed within a certain radius, as a function 
of that radius. Such panels summarize what we find in 
our simulation. As Fig. [T] indicates, we do indeed see an 
initial increase of angular momentum, but only between 
~60 km and ~250 km, and quasi-oscillatory behavior 
with a period near ~30 ms. The innermost zones that 
might be identified with the proto-neutron star are not 
spun up to periods faster than ~5— 10 seconds. At later 
times (the last 0.17 s after bounce), we note that the 
angular momentum that originally resided in this inner 
region is ejected with the exploding mass and follows the 
outward progress of the reenergized shock. This ejected 
angular momentum "bump" can also be seen clearly in 
Fig.|2j where we show the spatial profile of all three com- 
ponents of angular momentum for four different times 
near the end of our simulation. We address the subse- 



quent fate of this ejected angular momentum in [2.3 and 
^3j At the end of our simulation, very little angular mo- 
mentum is left in the central region of our grid, indicating 
that the inner core is left rotating very slowly. At the fi- 
nal timestep, the spin periods in the inner ten kilometers 
are all greater than ~10 seconds. Note that the flatness 
of the regions in Fig. [I] at large distances demonstrates 
that total angular momentum is very well conserved dur- 
ing the simulation. In CASTRO, angular momentum is 
not conserved by construction, so that its global conser- 
vation during this simulation is encouraging. Specifically, 
before and after bounce the total angular momentum of 
matter on the grid is conserved to better than ~1%, stay- 
ing very near "zero" for the duration of the calculation. 

Next, we divide our grid into thin spherical shells that 
cover the innermost 2000 km of the computational do- 
main. For the inner 250 km the width of the shells was 
chosen to be 1 km and for radii between 250 and 2000 km 
the width was increased to 10 km. For each of the shells 
we calculated the average specific angular momentum 
not the total enclosed angular momentum), for all three 
components of the angular momentum. Figure [3] depicts 
their associated spatial and temporal profiles. As dis- 
cussed above and can be seen in Figs, [I] and [2] regions 
interior to ~200 km experience some growth of specific 
angular momentum during the earlier stages after core 
bounce, whereas the specific angular momentum at larger 
radii remains close to zero until after the onset of explo- 
sion. As Fig. [3] makes clear, the sign of the specific angu- 
lar momentum for each component varies with time and 
radius, conserving total angular momentum after bounce 
to high precision. After the onset of explosion, angular 
momentum is seen to propagate outward with the ejected 
mass and the shock wave. 

An interesting feature is observed for radii between 60 
and 250 km. As time progresses in the simulation, the 
magnitudes of the specific angular momenta on shells in 
this region grow with time, indicating that in this region 
there is indeed some spin up by the spiral mode identified 
in BM07, Blondin & Shaw (2007), and Fernandez (2010). 



We proceed in [2.2 with a discussion of this phenomenon 
and the associated spin periods. 

2.2. Rotational Period and Frequency Profiles 

Having calculated the angular momentum contained 
in spherical shells at various radii, we proceed to extract 
the rotational profiles of those shells. Figure [4] shows 
the temporal evolution of the average spin frequency (/^, 
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in Hz) of each shell, for rotation around the x, y and z 
axes (the three panels in Fig. H). We see that shells with 
radii R < 60 km and R > 250 Km do not develop any sig- 
nificant rotational motion during the entire simulation. 
However, at the very early times after core bounce and 
for radii between ^60 km and ~250 km in the shocked, 
turbulent regions, some rotational motion does develop, 
which with time is amplified. By the end of the simu- 
lation, those shells have acquired spin frequencies that 
reach ~150 Hz at a radius of ~100 km, corresponding 
to a spin period of ~7 ms. Furthermore, we notice that 
although at early post-bounce times the sign of the spin 
vector of those shells kept flipping (as indicated by the 
positive and negative frequency values shown in Fig. [4|, 
at the end of the simulation neighboring shells tendecfto 
align their rotational motion. The left panel of Fig. [5] 
depicts the velocity field on a ball at radius 90 km near 
the end of the simulation. As this figure clearly shows, 
coherent fields of rotation have emerged in this region 
(see also the panels in Fig. [3. Near the same final time, 
the shock has been launchea into explosion and has the 
multi-turbule morphology portrayed in the right panel of 
Fig. [5] as an isoentropy surface. 

In Fig. [6] we show the temporal evolution of the magni- 
tude of the specific angular momentum (top) for various 
shells at different radii, and the corresponding magnitude 
of the rotational frequency (bottom) of those shells. As 
in Fig. |4j the development of rotational motion from ~60 
to ~250 km is obvious and the high spin rate of ^150 
Hz in this narrow range of radii is manifest. However, 
there is very little mass (~0. 01— 0.03 M©) and, hence, 
total angular momentum in this fastest spinning region, 
and this has major consequences for the spin evolution 
of the entire residual neutron star. 

2.3. Rotational Periods for the Final State 

Our simulation shows not only that the inner core isn't 
spun up significantly, but that the mass that does expe- 
rience an interesting degree of induced rotation is likely 
to be ejected with the subsequent blast. An interest- 
ing question arises associated with this ejected angular 
momentum "bump," shown in Figs, [I] and [2] If that 
angular momentum were to be accreted later on onto 
the central regions (in a "fallback" scenario for angular 
momentum), what would that imply for the terminal ro- 
tation rate of the nascent neutron star? We address this 
issue by applying various mass cuts to the final state of 
our simulation. We calculate the total angular momen- 
tum enclosed within a given spherical mass cut (given in 
column 2 of Table [T]) and then divide it by a represen- 
tative moment of inertia for the final-state neutron star 
(taken to be 2 x 10 45 g cm 2 and assuming solid-body 
rotation) to derive a final angular frequency. Our results 
are summarized in Fig. [7] and Table p] For a range of 
mass cuts between 1.2 and 1.6 M (baryonic) for our 
"1.9 x 10 52 ergs s _1 " explosion model the final inferred 
rotational periods vary between 1.2 and 47 seconds (as 
can be seen in Fig. [7]). The fastest inferred period of 
1.2 s is seen for a spherical mass cut of 1.53M0 baryonic. 
Note that many measured neutron star masses are near 
1.35 M gravitational, which translates into ~1.5 M 
baryonic. 

Apart for the magnitude of such final-state rotational 
periods, we would also like to know the direction of ro- 



tation and how that varies with the different spherical 
mass cuts. We define two angles to describe the direc- 
tion of the total angular momentum enclosed within a 
certain mass cut. Angle (0° ^ ^ 180°) is the polar 
angle, measured with respect to the y-axis, and angle (j) 
(0° ^ (j) ^ 360°) is the azimuthal angle, defined with 
respect to the z-axis. The two last columns of Table [I] 
show the values of and <p for the different mass cuts. 
The direction of the angular momentum varies little for 
spherical mass cuts from ~1.2 to ~1.4 M ). However, 
as Table [I] indicates, it does shows significant variation 
for mass cuts in the range 1.45 to 1.6 M . 

3. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the temporal and spa- 
tial evolution of the rotational profile of a newly-formed 
neutron star. The goal was to test the possibility that 
the remnant neutron star created in the context of core- 
collapse supernova could be induced to rotate rapidly 
due to the growth of a spiral, non-axisymmetric mode 
of the so-called SASI, even if the progenitor was initially 
non-rotating. We employed the 3D hydrodynamics code 
CASTRO to carry out the simulation, starting from the 
infall phase of a 15-M non-rotating progenitor, and fol- 
lowed the core bounce, convection, and explosion stages. 
The core, whose induced rotation was being studied, was 
included fully on the computational grid, a realistic EOS 
was used, and a simple neutrino heating and cooling al- 
gorithm enabled us to simulate in approximate fashion 
the turbulent post-bounce delay and subsequent explo- 
sive phases. Our main findings are: 

• We did not find that a significant spin rate was 
induced in the nascent neutron star. By the end 
of our simulation the inner core is seen to rotate 
very slowly, with spin periods no faster than ~5— 10 
seconds. 

• However, in the unstable gain region between the 
shock and the inner core (60 km < R < 250 km), 
counter-rotating shells can acquire enough specific 
angular momentum during the stall phase before 
explosion that they can achieve interestingly high 
spin rates, reaching frequencies around ~150 Hz 
(Fig. pi) and spin periods between ~6 and ~100 
ms. However, the mass in this region, and the cor- 
responding angular momentum are rather low. 

• At the later stages in the evolution of our simu- 
lation, the explosion and ejecta transport angular 
momentum outward behind the exploding shock 
wave. Even if some of this angular momentum were 
to fall back at later times, the minimum spin pe- 
riod that the residual neutron star could achieve is 
likely to be quite modest (Fig. [7] and Table [T]). 

Blondin & Mezzacappa (2007) employed a gamma-law 
equation of state, did not follow the hydrodynamics of 
the inner core, did not include in their simulations the 
collapse or the explosion phases, and did not include any 
neutrino interactions. Yet, they concluded that even a 
non-rotating progenitor could in principle leave behind a 
rotating remnant. In this qualitative conclusion we agree 
with Blondin & Mezzacappa (2007). 
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However, we used a general nuclear equation of state, 
included the entire core all the way to the center in the 
hydrodynamic simulation, included in a crude fashion the 
effects of neutrino heating, cooling, and electron capture, 
and followed both the collapse and the explosion phases. 
What we found was that the degree of induced central ro- 
tation of the residual proto-neutron star, after explosion 
and for any given degree of fallback, produces periods be- 
tween ~1 and ~10 seconds. Nevertheless, those authors 
suggest that a longer delay could lead to even greater net 
spin rates and see in their simulation a significant ramp 
up from ~0.4 to ~0.9 seconds after bounce. However, 
they concede that the spiral mode spin-up phenomenon 
should not continue beyond the onset of explosion. Since 
our explosion occurs near ^250 ms, one could argue that 
a longer delay could result in a larger induced spin. While 



this seems plausible, iWongwathanarat et al. (2010) have 



performed a new set ol approximate 3D simulations of 
core-collapse supernovae with non-rotating progenitors, 
focusing on the remnant neutron star kick velocities (see 
also Nordhaus et al. 2010b). They also show that after 
as much as ~1.4 seconds after bounce the induced spin 
rates are still fairly low, with final periods in the range 
~500 to ~1000 ms. Such long periods are approximately 
in accord with those we derive here in our "optimum" 
case. 

In summary, the increase in the spin rate from the op- 
timum derived from our simulation, which termina ted at 
^0.422 seconds, to that s een in the simulations of|Wong^| 



[wathanarat et al. (2010]), which terminated at ^1.4 sec- 
onds, was a factor ol ~2. The calculations performed by 
Blondin & Mezzacappa (2007) implied that the enhance- 
ment in delaying from ~0.422 seconds to ~0.9 seconds 
might be a factor of ~10 or more. We take this as a 
further indication that induced rotation by spiral modes 
may not be adequate to explain the observed (or inferred) 
rotation rates of most pulsars. Nevertheless, induced 
spin remains an intriguing, if sub-dominant, potential 
contribution to the overall angular momentum budget 
of nascent neutron stars. Where final neutron-star spin 
periods of seconds are observed, the spiral-mode spin-up 
mechanism may remain viable. 

One could argue that the delay to explosion we wit- 
ness is already rather long. We suggest that in order 
for the neutrino mechanism to be robust and to generate 
explosion energies sufficient to overcome the binding en- 
ergies of the progenitor mantle, while at the same time 
yielding explosion energies at infinity of ~10 51 ergs, the 
onset of explosion should probably be rather early after 
bounce. However, this has yet to be demonstrated. We 
emphasize that our models start with zero angular mo- 
mentum. A longer delay and total angular momentum 



conservation would naturally make it less and less likely 
that any significant net angular momentum would be left 
behind. This effect is implicit in Fig. [7| for which the 
residual net spin period jumps up on the right-hand side 
at greater mass cuts. Individual zones could initially be 
left with interesting spins, but others would have to be in 
the opposite directions required by net angular momen- 
tum conservation. At very late times, the residue would 
have to be in solid-body rotation at low net angular mo- 
mentum. 

In this paper, we have not explored various progenitor 
masses, nor a range of initial spin structures. Our goal 
was to determine whether a sizable rotation rate could 
be induced in the residual proto-neutron star when the 
progenitor itself was non-rotating. The behavior of the 
modest spiral mode we do see, and the overall role of 
rotation in the supernova phenomenon itself, remain to 
be fully mapped out. As computational capabilities im- 
prove, and codes acquire more physics and sophistica- 
tion, all these issues can and will be readdressed. 
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TABLE 1 

Results at the Final Stage Of The Simulation^) 
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/ 
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[xl0 46 g crn^- 1 ] 
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[Hz] 
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[°] 


1.20 


0.501 


2.50 


0.40 


97.6 


67.2 


1.25 


0.531 


2.37 


0.42 


92.1 


67.5 


1.30 


0.479 


2.83 


0.35 


89.7 


70.5 


1.35 


0.380 


3.30 


0.30 


96.7 


84.2 


1.40 


0.252 


4.99 


0.20 


111.6 


272.3 


1.45 


0.125 


10.17 


0.1 


152.5 


0.25 


1.50 


0.619 


2.02 


0.49 


128.6 


51.6 


1.55 


0.790 


1.59 


0.63 


21.6 


5.3 


1.60 


0.027 


46.89 


0.02 


125.5 


191.9 



( a ) At t = 0.422 s after core bounce; In column 1, we list various spherical mass cuts. The total angular momentum enclosed within that 
mass is given in column 2. If this enclosed angular momentum were to be left with the nascent neutron star, which we assume to have 
a moment of inertia of I ~ 2 x 10 45 g cm 2 , then the resulting rotational periods and frequencies would be those presented in columns 
3 and 4, respectively. The direction of the total enclosed angular momentum L for each case is represented by two angles given in the 
last two columns. The polar angle 6 (0° < 6 < 180°) is defined with respect to the y-axis and the azimuthal angle 4> (0° ^ <f> ^ 360°) is 
measured with respect to the z-direction. 
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Fig. 1. — Three-dimensional surface representation of the temporal evolution (after bounce) of the components (Li) of the total angular 
momentum (vertical axes) enclosed within a spherical radius R (in km). The three panels, from top to bottom, correspond to the x-, y-, and 
z-components of the enclosed angular momentum. The two-dimensional plot under each surface shows the projection of the upper surface 
onto the plane, with the color variation indicating the value of the projected point, going from black (negative) to white (positive) through 
a scale of red shades. We notice the initial increase of the angular momentum components in the very inner regions of our computational 
domain, and the ejection (with the ejected mass) of the bulk of angular momentum outwards during the later stages of the simulation. See 
text in §2.1| for further discussion. 
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Fig. 2. — Total angular momentum enclosed within a spherical radius R and its evolution with time after bounce, for four different 
timesteps. The three panels correspond to the three components of the angular momentum vector (x, y and z components from top to 
bottom, respectively). Lines of different color indicate different post-bounce times during the simulation. We notice again (as in Fig. fTj an 
initial increase of the angular momentum components between ~60 and ~200 km. With time, much of the angular momentum is ejected 
with the exploding mass. This can be seen here in the propagating angular momentum "bump." As discussed in the text and also in the 
captions to Tableland Fig. [7] this ejected angular momentum, if actually all accreted onto the neutron star at later time, would induce a 
rotational period no faster than 1.2 s. 
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Fig. 3. — Temporal evolution (in seconds after bounce) of the average specific angular momentum vertical axes) of spherical shells 
at various enclosed spherical radii (R, in centimeters). R is logaritmically distributed. From top to bottom, the three panels correspond 
to the x-, y-, and z-component of the specific angular momentum. To determine the averages, the computational grid was divided into 
thin spherical shells of width 1 km within the inner 250 km of the grid and of width 10 km for radii between 250 and 2000 km. The two- 
dimensional surfaces under each upper surface are projections of the upper surfaces, with the color variation indicating the value of the 
projected point, going from black (negative) to white (positive) through a scale of r ed s hades. The propagation outward at later times of 
the induced angular momentum is manifest in these panels as well. See the text in §2.1| for a discussion. 
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Fig. 4. — The temporal evolution of the spin frequency profiles of thin spherical shells at various spherical radii R (in centimeters, shown 
logarithmically). For each of the individual shells we have calculated the corresponding average rotational frequency, shown here on the 
vertical axis of each plot. Time is the post-bounce time in seconds. The three panels correspond to rotation around the x, y and z axes 
(from top to bottom). The positive and negative values of the frequency indicate opposite rotation directions. The computational grid was 
divided into thin spherical shells of width 1 km within the inner 250 km of the grid and of width 10 km for radii between 250 and 2000 km. 
We notice that for radii in the range R < 60 km and R > 250 km, the resulting frequencies remain near "zero" throughout the entire 
simulation. Interestingly, in the region 60 < R < 250 km, some rotational motion develops even during the early stages of the simula tion 
and amplifies during the very late stages, producing shell rotational frequencies that reach ~150Hz at a radius of 100 km. See text in §2.2| 
for a discussion. 
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Fig. 5. — Left: Velocity vector field profile on a spherical shell of radius 90 km towards the end of the simulation (~370 ms after core 
bounce). Both the size and color of the vectors represent the magnitude of the velocity field, with darker (lighter) shades corresponding 
to higher (lower) values. Right: 3D rendering of an isentropic surface with entropy equal to 10 kt baryon -1 . The colormap represents 
the magnitude of the entropy gradient on the surface, with values increasing going from orange to purple. The surface spans 1600 km and 
corresponds to a time 380 ms after core bounce. 
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Fig. 6. — Temporal evolution (in seconds) of the magnitude of the specific angular momentum (top) and of the magnitude of the rotational 
frequency (bottom) of spherical shells at various spherical radii R (in centimeters, shown logarithmically). As discussed in the caption to 
Fig. 4 and in the text, individual shells at radii 60 < R < 250 km develop significant rotational frequencies over time. However, outside 
this range of radii the rotation induced by non-axisymmetric instabilities is quite modest at all times during the simulation. 
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Fig. 7. — Rotation periods of inferred final-state neutron stars as a function of various possible mass cuts (in Mq), as described in the 
caption to Table 1 and in the text. To derive these approximate final spin periods, the enclosed total angular momentum is divided by a 
final neutron star moment of inertia of 2 x 10 45 g cm 2 . Solid-body rotation for the final state is assumed. The fastest induced spin up is 
seen to be for a mass cut of 1.53 Mq (baryonic mass), resulting in a slow rotational period of ~1.2 seconds. 



